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Abstract 

Spatial birth-and-death processes with a finite number of particles are obtained as unique 
solutions to certain stochastic equations. Conditions are given for existence and uniqueness 
of such solutions, as well as for continuous dependence on the initial conditions. The possibil¬ 
ity of an explosion and connection with the heuristic generator of the process are discussed. 
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Introduction 

This article deals with spatial birth-and-death processes which may describe stochastic dynamics 
of spatial population. Specifically, at each moment of time the population is represented as a 
collection of motionless points in We interpret the points as particles, or individuals. Existing 
particles may die and new particles may appear. Each particle is characterized by its location. 

The state space of a spatial birth-and-death Markov process on with finite number of 
points is the space of finite configurations over 

ro(M'') = {7]CR^: M < oo}, 

where \ri\ is the number of points of r]. 

Denote by the Borel cj-algebra on The evolution of a spatial birth-and-death 

process in admits the following description. Two functions characterize the development 
in time, the birth rate coefficient b : x ro(K‘^) —)■ [0;oo) and the death rate coefficient 

d : X ro(K‘^) —)■ [0;oo). If the system is in state rj G ro(K'^) at time t, then the probability 

that a new particle appears (a “birth”) in a bounded set B G ^{R'^) over time interval [t] t + At] 
is 

At J b{x,ri)dx + o{At), 

B 

the probability that a particle x £ r] is deleted from the configuration (a “death”) over time 
interval [t]t + At] is 

d{x, ri)At -|- o{At), 

and no two events happen simultaneously. By an event we mean a birth or a death. Using a 
slightly different terminology, we can say that the rate at which a birth occurs in B is b{x, ri)dx, 
the rate at which a particle x £ rj dies is d{x,ri), and no two events happen at the same time. 

Such processes, in which the birth and death rates depend on the spatial structure of the 
system as opposed to classical Z_|_-valued birth-and-death processes (see e.g. [KM59], [CG], 
[Har63, Page 116], [AN72, Page 109], and references therein), were first studied by Preston in 
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[Pre75]. A heuristic description similar to that above appeared already there. Our description 
resembles the one in [GK06]. 

The (heuristic) generator of a spatial birth-and-death process should be of the form 

LF{rj) = / h{x,rj)[F{rjiJ x) — F{rj)]dx + E d{x,r]){F{r]\x) - F{rj)), (1) 

for F in an appropriate domain, where rjU x and r]\x are shorthands for rj U {x} and rj \ {x}, 
respectively. 

Spatial point processes have been used in statistics for simulation purposes, see e.g. [MS94], 
[MW04, chapter 11] and references therein. For application of spatial and stochastic models in 
biology see e.g. [Lev03], [FOK“*'14], and references therein. 

To construct a spatial birth-and-death process with given birth and death rate coefficients, 
we consider in Section 2 stochastic equations with Poisson type noise 


r]t{B) = J I[o-b{x,v,-)]i^)dNi{x,s,u) 

Bx(0;t]x[0;oo] 

Zx (0;t] X [0;oo) 

where {rjt)t>o is a suitable ro(K'^)-valued cadlag stochastic process, the “solution” of the equation, 
I A is the indicator function of the set A, B G is a Borel set, A^i is a Poisson point processes 

on X M_|_ X M_|_ with intensity dx x ds x du, N2 is a Poisson point process on Z x M_|_ x M_|_ 
with intensity ^ x dr x dv, is the counting measure on Z'^, rjQ is a (random) initial finite 
configuration, 6, d ; x Fo(K'^) —)■ [0; 00 ) are functions that are measurable with respect to the 
product (T-algebra =^(M) x ,i^(ro(K)) and {xj} is some collection of points satisfying rjs C {xj} 
for every moment of time s (the precise definition is given in Section 1.3.1). We require the 
processes Ni,N 2 ,rio to be independent of each other. Equation (2) is understood in the sense 
that the equality holds a.s. for all bounded B G =^(M‘^) and t > 0. 

Garcia and Kurtz studied in [GK06] equations similar to (2) for infinite systems. In the earlier 
work [Gar95] of Garcia another approach was used: birth-and-death processes were obtained as 
projections of Poisson point processes. A further development of the projection method appears 
in [GK08]. Fournier and Meleard in [FM04] considered a similar equation for the construction 
of the Bolker-Pacala-Dieckmann-Law process with finitely many particles. 

Holley and Stroock [HS78] constructed a spatial birth-and-death process as a Markov family 
of unique solutions to the corresponding martingale problem. For the most part, they consider 
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a process contained in a bounded volume, with bounded birth and death rate coefficients. They 
also proved the corresponding result for the nearest neighbor model in with an infinite number 
of particles. 

Kondratiev and Skorokhod [KS06] constructed a contact process in continuum, with the 
infinite number of particles. The contact process can be described as a spatial birth-and-death 
process with 

6(x,r/) = A^a(x - y), (i(x,r/) = l, 

y&v 

where A > 0 and 0 < a G Under some additional assumptions, they showed existence 

of the process for a broad class of initial conditions. Furthermore, if the value of some energy 
functional on the initial condition is finite, then it stays finite at any point in time. 

In the aforementioned references as well as in the present work the evolution of the system in 
time via Markov process is described. An alternative approach consists in using the concept of 
statistical dynamics that substitutes the notion of a Markov stochastic process. This approach 
is based on considering evolutions of measures and their correlation functions. For details see 
e.g. [FKK12a], [FKK14], and references therein. 

There is an enormous amount of literature concerning interacting particle systems on lat¬ 
tices and related topics (e.g., [Lig85], [Lig04], [KL99], [Aldl3], [Fral4], [Spi77], etc.) Penrose in 
[Pen08] gives a general existence result for interacting particle systems on a lattice with local 
interactions and bounded jump rates (see also [Lig85, Chapter 9]). The spin space is allowed to 
be non-compact, which gives the opportunity to incorporate spatial birth-and-death processes 
in continuum. Unfortunately, the assumptions become rather restrictive when applied to con¬ 
tinuous space models. More specifically, the birth rate coefficient should be bounded, and for 
every bounded Borel set B the expression 


^ d(x,7?) 

should be bounded uniformly in r/, y G r(M'^). 

Let us briefly describe the contents of the article. 

In Section 1 we introduce give some general notions, definitions and results related to Markov 
processes in configuration spaces. We start with conhguration spaces, which are the state spaces 
for birth-and-death processes, then we introduce and discuss metrical and topological struc¬ 
tures thereof. Also, we present some facts and constructions from probability theory, such as 
integration with respect to a Poisson point process, or a sufficient condition for a functional 
transformation of a Markov chain to be a Markov chain again. 

In the second section we construct a spatial birth-and-death process (yt)t>o as a unique 
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solution to equation (2). We prove strong existence and pathwise uniqueness for (2). A key 
condition is that we require b to grow not faster than linearly in the sense that 

b{x,rj)dx < ci\ri\ + C 2 . (3) 

The equation is solved pathwisely, “from one jump to another”. Also, we prove uniqueness 
in law for equation (2) and the Markov property for the unique solution. Considering (2) with a 
(non-random) initial condition a G ro(M'^) and denoting corresponding solution by {ri{a,t))t>o, 
we see that a unique solution induces a Markov family of probability measures on the Skorokhod 
space [0; oo) (which can be regarded as the canonical space for a solution of (2)). 

When the birth and death rate coefficients b and d satisfy some continuity assumptions, the 
solution is expected to have continuous dependence on the initial condition, at least in some 
proper sense. Realization of this idea and precise formulations are given in Section 2.1. The 
proof is based on considering a coupling of two birth-and-death processes. 

The formal relation of a unique solution to (2) and operator L in (1) is given via the mar¬ 
tingale problem, in Section 2.2, and via some kind of a pointwise convergence, in Section 2.5. 

In Section 2.4 we formulate and prove a theorem about coupling of two birth-and-death pro¬ 
cesses. The idea to compare a spatial birth-and-death process with some “simpler” process goes 
back to Preston, [Pre75]. In [FM04] this technique was applied to the study of the probability 
of extinction. 

1 Configuration spaces and Markov processes: miscellaneous 

In this section we list some notions and facts we use in this work. 

1.1 Some notations and conventions 

Sometimes we write oo and -|-oo interchangeably, so that / ^ oo and / —)• -|-oo, or a < oo and 
a < -|-oo may have the same meaning. However, -|-oo is reserved for the real line only, whereas 
oo have wider range of applications, e.g. for a sequence {xnjneN C we may write Xn oo, 
n — 7- CO, which is equivalent to \xn\ -|-co. On the other hand, we do not assign any meaning 
to Xn —)■ -|-CO. 

In all probabilistic constructions we work on some probability space (H,^,P), sometimes 
equipped with a filtration of it- algebras. Elements of H are usually denoted as cj. 

The set is the complement of the set A C H: H \ A. We write [a; 6], [a; b) etc. for 

the intervals of real numbers. For example, (a; 6] = {x G M | a < x < 6}, — co < a < b < -|-co. 
The half line M+ includes 0: M+ = [0; co). 
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1.2 Configuration spaces 

In this section we introduce notions and facts about spaces of configurations, in particular, 
topological and metrical structures on r(M'^) as well as a characterization of compact sets of 
r(M‘^). We discuss configurations over Euclidean spaces only. 

1.1 Definition. For d G N and a measurable set A C the configuration space r(A) is defined 
as 

r(A) = {7 C A : I 7 n Ar| < +00 for any compact K C 

We recall that |A| denotes the number of elements of A. We also say that r(A) is the space 
of configurations over A. Note that 0 G r(A). 

Let Z_|_ be the set {0,1, 2,We say that a Radon measure /j, on (M'^, is a counting 

measure on if /r(A) G 1^+ for all A G When a counting measure n satisfies additionally 

t^{{x}) < 1 for all X G we call it a simple counting measure. 

As long as it does not lead to ambiguities, we identify a configuration with a simple counting 
Radon measures on as a measure, a configuration 7 G r(M'^) maps a set B ^ AS into | 7 ni?|. 
In other words, 7 = 

One equips r(M'^) with the vague topology, i.e., the weakest topology such that for all 
/ G C'c(K'^) (the set of continuous functions on with compact support) the map 

r(M'')9 7^(7,/) := J;/(x)gM 

a;G7 

is continuous. 

Equipped with this topology, r(M'^) is a Polish space, i.e., there exists a metric on r(M'^) 
compatible with the vague topology and with respect to which r(M'^) is a complete separable 
metric space, see, e.g., [KK06], and references therein. We say that a metric is compatible with 
a given topology if the topology induced by the metric coincides with the given topology. 

For a bounded B <zW^ and 7 G r(M'^), we denote 5 ( 7 , B) = min{\x—y\ : x,y G 705 , x / y}. 
Let Br{x) denote the closed ball in of the radius r centered at x. 

A set is said to be relatively compact if its closure is compact. The following theorem gives 
a characterization of compact sets in r(M'^), cf. [KK06], [HS78]. 

1.2 Theorem. A set F C r(M'^) is relatively compact in the vague topology if and only if 

sup{ 7 ( 5 n( 0 )) + 5 “^( 7 , 5 n( 0 ))} < 00 (4) 

-yGF 

holds for all n G N. 
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Proof. Assume that (4) is satisfied for some F C r(M‘^). In metric spaces compactness 
is equivalent to sequential compactness, therefore it is sufficient to show that an arbitrary se¬ 
quence contains a convergent subsequence in r(M'^). To this end, consider an arbitrary sequence 
{7n}nGN C F. The supremum sup7n(i?i(0)) is finite, consequently, by the Banach-Alaoglu 

n 

theorem there exists a measure ai G C{Bi{0))* ^here (7(51(0))* is the dual space of C{Bi{0)) 

^ and a subsequence {7n^^} C {jn} such that 7n^^|_Bi(o) (7(51(0))*. Furthermore, one 

may see that ai G r(5i(0)) (it is particularly important here that sup{(5“^(7, 5i(0))} < oo ). 

■y&F 

Indeed, arguing by contradiction one may get that ai{A) G for all Borel sets A, and Lemma 
1.5 below ensures that oi is a simple counting measure. 

Similarly, from the sequence we may extract subsequence {jn ^} C {7n^^} in such a way 
that 7 n ' converges to some a 2 G r(52(0)). Continuing in the same way, we will find a sequence 
of sequences such that qi™^^ ^ Om G r(5m(0)) and {qi™'"*'^^} C Consider now 

the sequence {qi^^lneN- For any m, restrictions of its elements to 5^(0) converge to am in 
r(5m(0)). Therefore, q^”^ —)■ a in r(M'^), where a = (Jan. 

n 

Conversely, if (4) is not fulfilled for some no G N, then we can construct a sequence {qnjneN C 
F such that either the first summand in (4) tends to infinity: 

qj(5„o(0)) oo,j ^ oo 

in which case, of course, there is no convergent subsequence, or the second summand in (4) tends 
to infinity. In the latter case, a subsequence of the sequence {qn|B„j,(o)}nGN niay converge to a 
counting measure (when all q„ are considered as measures). However, the limit measure can not 
be a simple counting measure. Thus, the sequence {qnjneN C F does not contain a convergent 
subsequence in r(M'^). □ 

We denote by CS{T(R'^)) the space of all compact subsets of r(M'^). 

1.3 Proposition. The topological space r(M'^) is not a - compact. 

Proof. Let {iLmjneN be an arbitrary sequence from (75'(r(M'^)). We will show that (JiL^ ^ 

n 

r(M‘^). To each compact Km we may assign a sequence q^\q^\ ... of positive numbers such 
that 


sup {q(5„(0)) +S ^(q,5„(0))} < q^\ 

There exists a configuration whose intersection with 5„(0) contains at least q^'^ + 1 points, 
for each n G N. This configuration does not belong to any of the sets {iLmjmGN) hence it can 
not belong to the union |J Km- D 
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Remark. Since r(M'^) is a separable metrizable space, Proposition 1.3 implies that r(M'^) 
is not locally compact. 

For another description of all compact sets in r(M‘^) we will use the set C of all 

positive continuous functions cj) satisfying the following conditions: 

1 ) (/>(x) = (/)(y) whenever |a:| = |y|, x,?/ G 

2) lim|3,|^oo = 0. 

For (j) ^ ^ we denote 

^ := / y}. 

1.4 Proposition, (i) For all c > 0 and (/> G <h 


Kr:= 


7 : JJ '^<t>ix,yh{dx)'y{dy) ^ cj G C'S’(r(M‘^)); 


(ii) For all K G C'S'(r(M^)) there exist (/> G <1> such that 





^^{x,y)'y{dx)-/{dy)} ^ 1 . 


Proof, (i) Denote 6n 
For 7 G Kc we have 


min 4>(x) > 0. 
xes„(o) 


c ^ 


'^{x,yh{dx)-f{dy) 


B„{0)xB„{0) 


and 


> jj <i>{x)(i>{y)I{x / y}-f{dx)j{dy) ^ 6'^7(R„(0))(7(Rn(0)) - 1) 
S„( 0 )xB„( 0 ) 

c ^ JJ 'i>{x,y)'y{dx)'y{dy) ^ ^ ^ -Bn(O)). 

S„(0)xS„(0) 


Consequently, 


sup 7 (Rn( 0 )) ^ 6'n\/c+ 1, 
yeKc 


and 


sup 5 \j,Bn{0)) ^ 

Vn 


It remains to show that Kc is closed, in which case Theorem 1.2 will imply compactness of Kc- 
The space r(M'^) is metrizable, therefore sequential closedness will suffice. Take 'jk G Kc, 7fc — t 7 
in r(M'^), k — )• 00. For n G N, let Tjj G C'c(M'^ x M*^) be an increasing sequence of functions such 
that ^ T, '^nix,y) = 'I'(x,y) for x,y G satisfying |x|,|y| < n, |x — y| > For such a 



sequence we have t for all x,y ^ x ^ y. For each / G C'c(M'^ x M'^), the 


map 


7 ] {t] X T], 




f{x,y)y{dx)y{dy) 


is continuous in the vague topology. Thus for all n G N, {'jk x 'yk,^n) —^(jx 7 ,^n}- Conse¬ 
quently, (7 X 7, < c, n G N, and by Fatou’s Lemma 


(7 X 7, 


^(x,y)7(dx)j(dy) = 


liminf 'i/n(x, y)7(dx)7(dy) < liminf / / ^n{x,y)7{dx)7{dy) < c. 

^ ^ J J 

' R‘*xR‘i 

To prove (ii), for a given compact set K C F(M'^) and a given function G denote 


an{K) := sup{7(.B„(0)) + (5 ^(7,.B„(0))} 

■y&K 


and 


bni(p) := sup \(p{x)\. 

|a;|>n 


Theorem 1.2 implies an{K) < 00, and we can estimate 


'^{x,y)7{dx)7{dy) = 


[B. 


n+1 


(0)\B„(0))x(b n+1 (0)\S„(0)) 


:5n+l 


(0)\B„ (0)) X (b„+i (0)\S„ (0)) 


(t){x)(l){y) ^^ / y}7{dx)7{dy) ^ 

\x - y\ 


^ JJ bl{an + l)7{dx)7{dy) ^ -h 1)^. 

(Bn+l(0)\Bn(0)) X (B„+l(0)\i?„(0)) 

Taking a function (p £ ^ such that 

6 1 


36„((/))(an + 1)'* < 


TT^ (n + 1)2 ’ 


we get 


sup { / / 'i/{x,y)7{dx)7{dy)} ^ !.□ 
7Gii' J J 

R'ixR'* 
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1.2.1 The space of finite configurations 

For A C the space ro(A) is defined as 

ro(A) ;= {r/ C A : |r/| < oo}. 

We see that ro(A) is the collection of all finite subsets of A. We denote the space of n-point 
configurations as rg”'^(A): 


r|,")(A) :={r,Gro(A) I |r,| =n}, n G N, 

and r® (A) := {0}- Sometimes we will write Fq instead of ro(K'^). Recall that we occasionally 
write r]\x instead of 77 \ {x} , r]{J x instead of 77 U {x}. 

To define a topological structure on ro(K'^), we introduce the following surjections (see, e.g., 
[KK02] and references therein) 


sym : | | —)■ Fo(M'^) 

n=0 

sym{{xi, ...,Xn)) = {xi, ...,Xn}, 

where 


(5) 


(]^d)n ...,Xn) G (K"*)"- I Xj G M"*, j = 1, ...,n,Xi / Xj,i / j}, ( 6 ) 

and, by convention, = {0}- 

The map sym produces a one-to-one correspondence between tt > 1, and the 

quotient space where is the equivalence relation on 


(X1,...,X„) (7/i,...,7/„) 

when there exist a permutation a : {1, ..., 7 r} ^ {1, ...,77} such that 


(®(t(1) ) ^(7(n)) ijjlt • •• jVn) • 

We endow Fq"'^(M'^) with the topology induced by this one-to-one correspondence. Equiv¬ 
alently, a set A C is open iff sym~^{A) is open in The space C 

we consider, of course, with the relative, or subspace, topology. As far as F®(M'^) = { 0 } is 
concerned, we regard it as an open set. 

Having defined topological structures on Fq"'^(M'^), rr > 0, we endow Fo(K'^) with the topology 


10 



of disjoint union, 


OO 

ro(M'') = □ (7) 

n=0 

N , , 

In this topology, a set K C ro(K^) is compact iff K C |J Fq (M‘^) for some iV G N and for 

n=0 

each n < N the set K n is compact in rQ"'^(M‘^). A set Kn C rQ"'^(M‘^) is compact iff 

sym~^{Kn) is compact in We note that in order for Kn to be compact, the set sym~^Kn, 

regarded as a subset of should not have limit points on the diagonals, i.e. limit points 

from the set \ 

Let us introduce a metric compatible with the described topology on ro(M'^). We set 


dist{(, rj) 


1 A de^cKC,??), ICI = l??l, 

1 , otherwise. 


Here dEucliCid) the metric induced by the Euclidean metric and the map sym'. 


dEuci{C,d) = - y\ -x e sym ^C,,y ^ sym ^r/}, (8) 

where \x — y\ is the Euclidean distance between x and y, sym~^rj = sym~^{{'r]}). In many 
aspects, this metric resembles the Wasserstein type distance in [RS99]. The differences are, dist 
is bounded by 1 and it is defined on ro(M'^) only. 

Note that the metric dist satisfies equalities 


dist{(Ux,riUx) = dist{C,ri) (9) 

for C,!] ^ ro(M'^), x G x ^ C,ri, and 

dist{C\x,ri\x) = dist{(,r]), (10) 

X £ C,r]. We note that the space ro(M'^) equipped with this metric is not complete. Nevertheless, 
ro(K'^) is a Polish space, i.e., ro(K'^) is separable and there exists a metric p which induces the 
same topology as dist does and such that ro(K'^) equipped with p is a complete metric space. To 
prove this, we embed into the space of n-point multiple configurations, which 

we define as the space of all counting measures rj on with = n. Abusing notation, 

we may represent each rj G as a set {xi, ...,Xn}, where some points among Xj G 

may be equal (recall our convention on identifying a configuration with a measure; as a measure, 

n 

r] = 6x )■ One should keep in mind that {xi,..., x^} is not really a set here, since it is possible 
that Xi = Xj for i ^ j, i,j G {1, The representation allows us to extend sym to the map 
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sym : 


( 11 ) 


y {R^)^ ^ r^Q\R^) 

m=0 

sym{{xi, ...,Xn)) ■■= {xi, ...,Xn}, 

and define a metric on for C,r] G fQ"'^(M^) we set dist{C,ri) = l/\dEud{Ci'n): dEucliC^d) 

is the metric induced by the Euclidean metric and the map sym: 

dEudiC, f]) = inf{|x -y\ :x G sym~^C, V G (12) 

The metrics dist and dist coincide on rQ”^(M‘^) x rQ”^(M‘^) (as functions). Furthermore, one 
can see that (rg”^(M‘^), dist) is a complete separable metric space, and thus a Polish space. The 
next lemma describes convergence in (compare with Lemma 3.3 in [KK06]). 

1.5 Lemma. Assume that rf^ y in rg"'^(M'^), and let y = {xi, ...,x„}. Then y^, m G N, may 
be numbered, r/™ = {x™, ...,x™}, in such a way that 

X™ Xi, m —)■ oo 

in R'^. 

Proof. The inequality dist{y^,y'^) < e implies existence of a point from y^ in the ball 
B^{xi) of radius e centered at Xj, i G {1, ...,n}. Furthermore, in the case when x* is a multiple 
point, i.e., if Xj = Xi for some j ^ i, then there are at least as many points from y"^ in B^{xi) 
as y{{xi]). Observe that, for e < i inf{|x — y\ : y{{x}),y{{x}) > 1} A 1, we have in the previous 
sentence “exactly as many” instead of “at least as many”, because otherwise there would not be 
enough points in y^. The statement of the lemma follows by letting e ^ 0. 

1.6 Lemma. ro(K^) is a Polish space. 

Proof. Since ro(M'^) is a disjoint union of countably many spaces rQ"^^(M'^), it suffices to 
establish that each of them is a Polish space. To prove that rQ”^(M'^) is a Polish space, n G N, 
we will show that it is a countable intersection of open sets in a Polish space fQ’^^(M'^). Then we 
may apply Alexandrov’s theorem: any Gs subset of a Polish space is a Polish space, see §33, VI 
in [Kur66]. 

To do so, denote by the closed ball of radius m in with the center at the origin. 
Define Fm '■= {y G fg”^(M'’*) | y{{x}) > 2 for some x G B^} and note that 


12 



m=l 

Since fQ”'^(M'^) is Polish, it only remains to show that is closed in This is an 

immediate consequence of the previous lemma. 


1.2.2 Lebesgue-Poisson measures 

Here we dehne the Lebesgue-Poisson measure on ro(K'^), corresponding to a non-atomic Radon 
measure a on Our prime example for a will be the Lebesgue measure on For any 
n G N the product measure o'®"' can be considered by restriction as a measure on The 

projection of this measure on Pq”^ via sym we denote by so that 


= a®^{sym-^A), A G 

On P® the measure is given by a^^\{0}) = 1. The Lebesgue-Poisson measure on 
(ro(K^), .^^(Po(M'^))) is dehned as 


V:=E 


rO' 


in) 


n=0 


n\ 


(13) 


The measure is hnite iff a is hnite. We say that a is the intensity measure of Xo 


1.2.3 The Skorokhod space 

For a complete separable metric space {E, p) the space of ah cadlag Pi-valued functions 
equipped with the Skorokhod topology is a Polish space; for this statement and related dehni- 
tions, see, e.g., Theorem 5.6, Chapter 3 in [EK86]. Let p£) be a metric on De compatible with 
the Skorokhod topology and such that {De,Pd) is a complete separable metric space. Denote 
by (V{De), Pp) the metric space of probability measures on ^^{De), the Borel a - algebra of 
De, with the Prohorov metric, i.e. for P,Qg V{De) 


Pp(P, Q) = inf{e > 0 : P{F) < Q{F^) + e for ah F G ^{De)} (14) 

where 

F’^ = {x G De ■■ pd{x,F) < s}. 

Then {V{De), Pp) is separable and complete; see, e.g., [EK86], Section 1, Chapter 3, and 
Theorem 1.7, Chapter 3. The Borel c-algebra AS{De) coincides with the one generated by the 
coordinate mappings; see Theorem 7.1, Chapter 3 in [EK86]. In this work, we mostly consider 
Dr(j(]Rd)[0; T] and Dp(]gd)[0;T] endowed with the Skorokhod topology. 
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1.3 Integration with respect to Poisson point processes 

We give a short introduction to the theory of integration with respect to Poisson point processes. 
For construction of Poisson point processes with given intensity, see e.g. [Kal02, Chapter 12], 
[Kin93], [RY05, Chapter 12, § 1] or [IW81, Chapter 1, § 8,9]. All definitions, constructions and 
statements about integration given here may be found in [IW81, Chapter 2, § 3]. See also [GS79, 
Chapter 1] for the theory of integration with respect to an orthogonal martingale measure. 

On some filtered probability space (P, {^}t> 0 , P), consider a Poisson point process N 

on M_|_ X X X M_|_ with intensity measure dt x /3{dx) x du, where X = or X = We require 

the filtration {^}t>o to be increasing and right-continuous, and we assume that is complete 

under P. We interpret the argument from the first space M+ as time. For X = the intensity 
measure (3 will be the Lebesgue measure on for X = Z'^ we set jS = where 

#A = \A\, Ag'^(Z‘^). 

The Borel cj-algebra over Z'^ is the collection of all subsets of Z^, i.e. = 2^^*. Again, as 

is the case with configurations, for X = we treat a point process as a random collection of 
points as well as a random measure. 

We say that the process N is called compatible with > 0) if X is adapted, that is, all 

random variables of the type N{Ti, U), Tf G =^([0; t]), U G 1 ^(X x M+), are ^j-measurable, and 
all random variables of the type N(t + h,U) — N(t, U), h > 0,U £ =^(X x M+), are independent 
of X(f, 17) = X([0;f],17). For any U G x M_|_) with (/3 x 1){U) < oo, I is the Lebesgue 
measure on the process (X([0;t],17) — t(3 x l{U),t > 0) is a martingale (with respect to 
> 0); see [IW81, Lemma 3.1, Page 60]). 

1.7 Definition. A process / : M+ x X x M_|_ x 17 —)■ M is predictable, if it is measurable with 
respect to the smallest a - algebra generated by all g having the following properties: 

(i) for each f > 0, {x,u,lj) e-)- g(t, x,u,lj)) is =^(X x M+) x measurable; 

(ii) for each {x,u,u}), the map t e-)■ g{t,x,u,uj)) is left continuous. 

For a predictable process / G L^([0;T] x X x M_|_ x 17), t G [0;T] and U G 7^(X x M_|_) we 

define the integral It{f) = f f{s,x,u,uj)dN(s,x,u) as the Lebesgue-Stieltjes integral with 

[0;i]xC/ 

respect to the measure N: 


j fis,x,u,uj)dN{s,x,u)= ^ f{s,x,u,uj). 

[0;i]xt/ s<t,{s,x,u)GN 

This sum is well defined, since 
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E E \f{s,X,U,Uj)\ = J \f{s,x,u,u)\ds/3{dx)du < oo 

s<t,{s,x,u)eN 

We use dN{s,x,u) and N{ds,dx,du) interchangeably when we integrate over all variables. The 
process /*(/) is right-continuous as a function of t, and adapted. Moreover, the process 


Itif) 


J f {s, X, u, u})[dN(s, X, u) — dsl3{dx)du] 

[0-,t]xU 


is a martingale with respect to (^t,t > 0), [IW81, Page 62]. Thus, 


E J f{s,x,u,u))dN{s,x,u)=E J f{s,x,u,io)dsl3{dx)du. (15) 

[0;t]x(7 [0-,t]xU 

This equality will be used several times throughout this work. 


1.8 Remark. We can extend the collection of integrands, in particular, we can define 

f f(s,x,u,uj)dN(s,x,u) for / satisfying 
[0;t]xU 


E 


{\f{s,x,u,uj)\ A l)dsj3{dx)du < oo. 


[0;t]x(7 

However, we do not use such integrands. 


The Lebesgue-Stieltjes integral is defined w-wisely and it is a function of an integrand and 
an integrator. As a result, we have the following statement. The sign = means equality in 
distribution. 


1.9 Statement. Let he Poisson point processes defined on some, possibly different, proba¬ 
bility spaces, and let at be integrands, k = \,2, such that integrals f akdMk are well defined. If 
{ai, Ml) = {a 2 , M 2 ), then 


aidMi= / a2dM2. 


The proof is straightforward. 


1.3.1 An auxiliary construction 

Let be the counting measure on [0,1], i.e. 

#C = |C|, CG.^^([0;1]). 


The measure is not cj-finite. For a cadlag ro(K'^)-valued process {'nt)tG[ 0 -,oo]i adapted to 
{^i}te[ 0 ;oo]) we would like to define integrals of the form 


15 



J I{x€Bn'nr-}f{x,r,v,uj)dN 2 {x,r,v) (16) 

X [0;oo] X [0;oo) 

where B is a bounded Borel subset of / is a bounded predictable process and N 2 is a Poisson 
point process on x [0;T] x [0; 00) with intensity ^ x dr x dv, compatible with {^t}t<E[o-oo\- 
We can not hope to give a meaningful definition for an integral of the type (16), because of 
the measurability issues. For example, the map 

n —)■ M, 

w N2{u{uj), [ 0 ; 1 ], [ 0 ; 1 ]), 

where u is an independent of N 2 uniformly distributed on [0; 1] random variable, does not have 
to be a random variable. Even if it were a random variable, some undesirable phenomena would 
appear, see, e.g., [Pod09]. 

To avoid this difficulty, we employ another construction. A similar approach was used in 
[FM04]. If we could give meaningful definition to the integrals of the type (16), we would expect 

j A V, Lo)dN 2 {x, r, v)- 

R^x [0;t] X [0;oo) 

-^{xeBnr,,_}/(2;, A V, uj)#{dx)drdv 

R^x [0;t] X [0;oo) 

to be a martingale (under some conditions on / and B). 

Having this in mind, consider a Poisson point process A ^2 on Z x M_|_ x M_|_ with intensity 
^ X dr X dv, defined on (H, H) (here ^ denotes the counting measure on Z. This 

measure is cr-finite). We require N 2 to be compatible with {^}t>o- Lot (^t)te[ 0 ,oo] Le an adapted 
cadlag process in ro(M'^), satisfying the following condition: for any T < 00, 

Rt = \ IJ Vtl <00 a.s. (17) 

te[0;r] 

The set R^o ■= Ute [0 oo] ^t most countable, provided (17). Let ^ be the lexicographical 
order on We can label the points of r]Q, 

r]o = {X0,X-1, ...,X-q}, Xo 4 X-l 4 ... 4 X-g. 

There exists an a.s. unique representation 

Roo\'no = {xi,X2,...} 


16 



such that for any n,m G N, n < m, either inf{s : Xn G rjs} < inf{s : Xm G rjg}, or inf{s : Xn G 

s>0 s>0 s>0 

rjs} = inf{s ; Xm G rjs} and Xn =4 Xm- In other words, as time goes on, appearing points are 
s>0 

added to {xi,X 2 , •••} in the order in which they appear. If several points appear simultaneously, 
we add them in the lexicographical order. 

For the sake of convenience, we set x-i = A, i < —q — 1, where A ^ Z. We say that the 
sequence {...,x_i,xi, X 2 , •..} is related to (??t)te[ 0 ;oo]- 

For a predictable process / G x M+ x M+ x Q) and B G consider 


J I{xierir-nB}fixi,r,v,uj)dN2ii,r,v). (18) 

Zx (£i ;t 2 ] X [0;oo) 

Assume that is bounded for some T > 0. Then, for a bounded predictable / G x 

X M+ X Q) and B G the process 


I{xierjr-nB}f{xi, r, v, uj)dN 2 (i, r, v) 


Zx (0;^] X [0;oo) 


Zx {0;t] X [0;oo) 

is a martingale, cf. [IW81, Page 62]. 


J I{xi&vr-nB}fixi, r, V, u!)#{di)drdv 


1.3.2 The strong Markov property of a Poisson point process 

We will need the strong Markov property of a Poisson point process. To simplify notations, 
assume that A is a Poisson point process on M+ x with intensity measure dt x dx. Let N be 
compatible with a right-continuous complete filtration {^t}t>o, and r be a finite a.s. {4^t}t>o- 
stopping time (stopping time with respect to {<^t}t>o)- Introduce another Point process N on 
M+ X M"*, 


A([0; s] X t/) = A((r; t + s]xU), U G ^(M'^). 

1.10 Proposition. The process N is a Poisson point process with intensity dt x dx, independent 

of 

Proof. To prove the proposition, it is enough to show that 

(i) for any b > a > 0 and open bounded U C N{{a; b), U) is a Poisson random variable 
with mean (b — a)/3{U), and 

(ii) for any bk > Ok > 0, k = l,...,m, and any open bounded Uk C M^, such that {{afjbi) x 
Ui)r\{{aj;bj)xUj) = 0,i /j,the collection {N{{ak;bk)xUk)}k=i,m is a sequence of independent 
random variables, independent of 
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Indeed, N is determined completely by values on sets of type (b — a)l3(U), a, b, U as in (i), 
therefore it must be an independent of Poisson point process if (i) and (ii) hold. 

Let Tn be the sequence of {^j}j>o-stopping times, ^ on {r G ^]}, A: G N. Then 
Tn -I T and Tn — T < The stopping times Tn take only countably many values. The process 
N satisfies the strong Markov property for Tn- the processes Nn, defined by 

Afn([0; s] X U) := N{{Tn] Tn + s] X U), 

are Poisson point processes, independent of ^t„- To prove this, take k with P{Tn = ^} > 0 
and note that on {r^ = ^}, Nn coincides with process the Poisson point process N_k_ given by 

%([0; 5] X [/) := iv(( A; A + s] X I/)) , [7 G 

Conditionally on {r^ = ^}, N_k_ is again a Poisson point process, with the same intensity. 

Furthermore, conditionally on {Tn = -M;}, N k is independent of ^ fc , hence it is independent 

^ 2 " 2 " 

of c ,^_k_ . 

2 " _ _ 

To prove (i), note that Nn{{a;b) x 17) —)• N({a]b) x U) a.s. and all random variables 
Nn{{a] b) X U) have the same distribution, therefore N{{a; b) x [/) is a Poisson random variable 
with mean (b — a)X(U). The random variables Nn{{a-,b) x U) are independent of hence 
N{{a;b) x U) is independent of too. Similarly, (ii) follows. □ 

Analogously, the strong Markov property for a Poisson point process on M+ x N with intensity 
dt X ^ may be formulated and proven. 

1.11 Remark. We assumed in Proposition 1.10 that the filtration compatible with 

N, is right-continuous and complete. To be able to apply Proposition 1.10, we should show that 
such nitrations exist. 

Introduce the natural filtration of N, 

^0 = a{NkiC,B),B£ £ ^i[0;t])}, 

and let be the completion of under P. Then N is compatible with We claim 

that {-^t}t>o, defined in such a way, is right-continuous (this may be regarded as an analog of 
Blumenthal 0—1 law). Indeed, as in the proof of Proposition 1.10, one may check that Na is 
independent of =^a+- Since = (T{Na) V (T{Na) and are independent and ^a+ C ^oo; 
one sees that ^a-i- C ^a- Thus, ,^a+ = ^a- 

1.12 Remark. We prefer to work with right-continuous complete filtrations, because we want 
to ensure that there is no problem with conditional probabilities, and that the hitting times we 
will consider are stopping times. 
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1.4 Miscellaneous 


When we write ^ ~ Exp{X), we mean that the random variable ^ is exponentially distributed 
with parameter A. 

1.13 Lemma. If a and /3 are exponentially distributed random variables with parameters a and 
b respectively (notation: a ~ Exp{a), /3 ~ Exp{b) ) and they are independent, then 


Indeed, 


P{a < 13} 


a 

a + b 


fOO I'OO 

P{a <(3}= aP{x < /3}e-“^ = a / 

Jo Jo a + b 

Here are few other properties of exponential distributions. If ^i, ^ 2 ; •••) are independent ex¬ 
ponentially distributed random variables with parameters ci,...,Cn respectively, then min 
is exponentially distributed with parameter ci c^- Again, the proof may be done by direct 

computation. If are independent exponentially distributed random variables with pa¬ 

rameter c and ai,Q; 2 ) ••• is an independent sequence of independent Bernoulli random variables 
with parameter p G (0; 1), then the random variable 

e 

^ ^ 0 = mm{k G N : Ofc = 1} 

i=l 

is exponentially distributed with parameter The random variable f is the time of the first 
jump of a thinned Poisson point process with intensity c. The statement about the distribution 
of ^ is a consequence of the property that the independent thinning of a Poisson point process 
with intensity A is a Poisson point process with intensity pX, see [Kal02, Theorem 12.2,(iv)]. 

We will also need the result about finiteness of the expectation of the Yule process. A Yule 
process (Zt)t>o is a pure birth Markov process in Z_|_ with birth rate pn, p > 0, n £ Z_|_. That 
is, if Zt = n, then a birth occur at rate pn, i.e. 

P{Zt+At - Zt = I \ Zt = n} = pn + o{At). 

For more details about Yule processes see e.g. [AN72, Chapter 3], [Har63, Chapter 5], [ArnOG] 
and references therein. Let {Zt{n))t>o he a Yule process started at n. The process {Zt{n))t>o 
can be considered as a sum of n independent Yule processes started from 1, see e.g. [ArnOG]. 
The expectation of Zt{l) is finite and EZt{l) = see e.g. [AN72, Chapter 3, Section 6] or 
[Har63, Chapter 5, Sections G,7]. Consequently, if {Zt)t>o is a Yule process with EZq < oo, then 
EZt < oo and EZt = EZoe^^. 
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Here are some other properties of Poisson point processes which are used throughout in the 
article. If is a Poisson point process on M_|_ x x M_|_ with intensity ds x dx x du, then a.s. 

Vx G M"* : N{R+ x {x} x M+) < 1. (19) 

Put differently, no plane of the form M+ x {x} x M+ contains more than 1 point of N. Using 
the (T-additivity of the probability measure, one can deduce (19) from 

Vx G : A1([0;1] X {x} X [0;1]) < 1. (20) 

We can write 

|vx G M"* ; iV([0; 1] x {x} x [0; 1]) < 1 

D jvfe G {0,1, ...,n - 1} : iV([0; 1] x [^; x [0; 1]) < l|, 

and then we can compute 

pjvfc € {0,1. ...,n - 1) : JV([0i 1| x |^; X [0; 1]) < l| 

= (P{W(|0; 1] X [0; i] X |0; 1]) <!})”= (exp(-i)|l + i|)" = (l - = 1 - o(i). 

Thus, (20) holds. 

Let G L^(M‘^), 'ip >0. Consider the time until the first arrival 

r = inf{t > 0 : J I[o;^(j;)](M)iV((is, dx, drt) > 0}. (21) 

[0;t]xR‘*xR+ 

The random variable r is distributed exponentially with the parameter HV’llii- From (19) we 
know that a.s. 


N{{t} X X M+) = lV({(r,x,u) | x G G [0;V'(a:)]}) = 1 
Let Xr be the unique element of defined by 


Then 


A^({r} X {xt} X M+) = 1. 


P{x.^ G B} 


JB'ipix)dx 

V'(x)dx’ 


B G ^(M'^). 


( 22 ) 
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1.5 Pure jump type Markov processes 

In this section we give a very concise treatment of pure jump type Markov processes. Most of 
the definitions and facts given here can be found in [Kal02, Chapter 12]; see also, e.g., [GS75, 
Chapter 3, § Ij. 

We say that a process X = {Xt)t>o in some measurable space {S,S) is of pure jump type if 
its paths are a.s. right-continuous and constant apart from isolated jumps. In that case we may 
denote the jump times of X by ri,r 2 ,..., with understanding that = oo if there are fewer 
that n jumps. The times are stopping times with respect to the right-continuous filtration 
induced by X. For convenience we may choose X to be the identity mapping on the canonical 
path space When X is a Markov process, the distribution with initial 

state X is denoted by Px, and we note that the mapping x i—)■ Px{A) is measurable in x, X G Q. 

Theorem 12.14 [Kal02] (strong Markov property, Doob) A pure jump type Markov process 
satisfies strong Markov property at every stopping time. 

We say that a state x G S' is absorbing if Px{X = x} = 1. 

Lemma 12.16 [Kal02] If x is non-absorbing, then under Px the time ri until the first jump 
is exponentially distributed and independent of Orj^X. 

Here 9t is a shift, and O^X defines a new process, 

er,X{s)= X{s + Tl). 

For a non-absorbing state x, we may define the rate function c(x) and jump transition kernel 
fi{x,B) by 


c(x) = (ExTi) \ = Px{Xr^ G B}, X G S', B eS. 

In the sequel, c(x) will also be referred to as jump rate. The kernel C/U is called a rate kernel. 

The following theorem gives an explicit representation of the process in terms of a discrete¬ 
time Markov chain and a sequence of exponentially distributed random variables. This result 
shows in particular that the distribution Px is uniquely determined by the rate kernel c^. We 
assume existence of the required randomization variables (so that the underlying probability 
space is “rich enough”). 

Theorem 12.17 [Kal02] (embedded Markov chain) Let X be a pure jump type Markov process 
with rate kernel cp.. Then there exists a Markov process Y on Z_|_ with transition kernel jj, and 
an independent sequence of i.i.d., exponentially distributed random variables 71 , 72 ,... with mean 
1 such that a.s. 


Xt — Yji, t G [tji, Tfi-i-i), n G Z^, 


(23) 
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where 


n G Z+. 


(24) 


Tn 


E 


k=l 


Ik 

c(n-i)’ 


In particular, the differences between the moments of jumps Tn+i — Tn of a pure jump type 
Markov process are exponentially distributed given the embedded chain Y, with parameter 
ciYn). If c(Y’fc) = 0 for some (random) k, we set = oo for n > A: +1, while are not defined, 
n > A: + 1. 

Theorem 12.18 [Kal02] (synthesis) For any rate kernel on S with fi{x, {x}) = 0, consider 
a Markov chain Y with transition kernel fi and a sequence 71 , 72 , ••• of independent exponentially 
distributed random variables with mean 1, independent of Y. Assume that cjFT ~ 
for every initial distribution for Y. Then (23) and (24) define a pure jump type Markov process 
with rate kernel cfi. 

Next proposition gives a convenient criterion for non-explosion. 

Proposition 12.19 [Kal02] (explosion) For any rate kernel cp. and initial state x, let (Yn) 
and (xn) be such as in Theorem 12.17. Then a.s. 


r„ —)• 00 





= 00 . 


In particular, —)• 00 a.s. when x is recurrent for (Yn). 


(25) 


1.6 Markovian functions of a Markov chain 

Let {S, SS{S)) be a Polish (state) space. Consider a (homogeneous) Markov chain on {S, FS{S)) as 
a family of probability measures on S^. Namely, on the measurable space (fl, = {S^ ,FS{S^)) 
consider a family of probability measures such that for the coordinate mappings 


Xn-.^^ S, 

^n(si, S2,'") 

the process X = {Xn}nGZ+ is a Markov chain, and for all s G S' 

= s} = 1, 

Ps{Xn+mj G Aj, j = I Mn) — Px^ {Xmj G Aj, j = 1, ..., /ci}. 

Here Aj G mj G N, Ai G N, = o'lXi,..., Xn}. The space S is separable, hence 

there exists a transition probability kernel Q : S x AS{S) —)■ [0; 1] such that 
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Q{s, A) = Ps{Xi £ A}, s £ S, Ag Ag{S). 


Consider a transformation of the chain X, = f{Xn), where / : 5 —)■ Z_|_ is a Borel- 
measurable function, with convention =^(Z_|_) = 2^+. In the future we will need to know when 
the process Y = {Yn}z+ is a Markov chain. A similar question appeared for the hrst time in 
[BR58]. 

A sufficient condition for Y to be a Markov chain is given in the next lemma. 

1.14 Lemma. Assume that for any bounded Borel function h : S ^ S 

Esh{Xi) = Eqh{Xi) whenever f{s) = f{q), (26) 

Then Y is a Markov chain. 

Remark. Condition (26) is the equality of distributions of Xi under two different measures, 
Pg and Pq. 

Proof. For the natural filtrations of the processes X and Y we have an inclusion 

D n G N, (27) 

since M is a function of X. For k £N and bounded Borel functions hj : Z+ -£■ W, j = 1,2, ...,k 
(any function on Z_|_ is a Borel function), 


Eg 


llhq{Y^+j)\^, 

i=i 


EX^UhqifiXq)) 

i=i 


/ Q{xo,dxi)hi{f{xi)) / Q{xi,dX 2 )h 2 {f{x 2 ))... / Q{Xn-l,dXn)hn{f{Xn)) 

Is Js Js 


(28) 


XQ — Xn 


To transform the last integral, we introduce a new kernel: for y £ f{S) chose x £ S with 
f{x) = y, ans then for B C Z_|_ define 


Qiy,B) = Qix,f-^iB)y, 


(29) 


The expression on the right-hand side does not depend on the choice of x because of (26). To 
make the kernel Q defined on Z_|_ x .i^(Z_|_), we set 


Q{y,B) = I{oeB}, y i f{S). 

Then from the change of variables formula for the Lebesgue integral it follows that the last 


23 



integral in (28) allows the representation 


I / Q(^f (^Xn—i), dZfi^hn^Zfi). 

Js Jz+ 

Likewise, we set Zn-i = f{xn-i) in the next to last integral: 

I Q(^Xfi—2t dXn—l)hn(^f (^Xfi—l)) I Q{Xn—l, dXn')hii{f {Xn)) — 

Js Js 

I (5(Xyj_2, l)^n(/(^n—l)) / Q(/(^n—l); ^^'2'n)^n(^n) 

Js Jz+ 

I Qif iXn—2^ T dZn—l^hn{Zn—l^ I Q{Zn—l, dZfi)hn{Zfi} ■ 

J J 

Further proceeding, we get 


/ Q{xo,dxi)hi{f{xi)) / Q{xi,dX 2 )h 2 {f{x 2 )):. / Q{Xn-l,dXn)hn{f{Xn)) = 

Js Js Js 

/ Q{zo,dzi)hi{zi) / Q{zi,dZ 2 )h 2 {z 2 )... / Q{Zn-l,dZn)hn{Zn), 

■J ^_l_ J ^_l_ J 

where zo = f{xo). 

Thus, 






i=i 


IZ+ 


Q{f{Xo),dzi)hi{zi) / Q{zi,dZ2)h2{z2)... / Q{Zn-l,dZn)hn{Zn 


This equality and (27) imply that y is a Markov chain. 

1.15 Remark. The kernel Q and the chain f{Xn) are related: for all s G 5, n, m G N and 
M C N, 


Ps{f{Xn+i) G M I f{Xn) =m} = Q{m,M) 

whenever Ps{f {Xn+i) = m} > 0. Informally, one may say that Q is the transition probability 
kernel for the chain {f(Xn)}n£Z+- 

1.16 Remark. Clearly, this result holds for a Markov chain which is not necessarily defined on 
a canonical state space, because the property of a process to be a Markov chain depends on its 
distribution only. 
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2 A birth-and-death process in the space of finite configurations: 
construction and basic properties 

We would like to construct a Markov process in the space of finite configurations ro(K'^), with 
a heuristic generator of the form 

LF{r]) = / b{x,r])[F{r]U x) — F{r])]dx+ '^^d{x,r]){F{r]\x) — F{r])). (30) 

for F in an appropriate domain. We call the functions 6 : x ro(K'^) —)■ [0;oo) and d : 

X ro(K'^) —)• [0; oo) the birth rate coefficient and the death rate coefficient, respectively. 
Theorem 2.16 summarizes the main results obtained in this section. 

To construct a spatial birth-and-death process, we consider the stochastic equation with 
Poisson noise 


i?x(0;t]x[0;oo] 

(31) 

I{xievr-nB}I[o-,d(xi,Vr-)] {v)dN-2{i, r, v) -h m{B), 

Zx (0;t] X [0;oo) 

where {r]t)t>o is a suitable cadlag ro(K‘^)-valued stochastic process, the “solution” of the equation, 
B G is a Borel set, Ni is a Poisson point process on x M_|_ x M+ with intensity 

dx X ds X du, N ‘2 is a Poisson point process on Z x M_|_ x M_|_ with intensity ^ x dr x dv ] rjo 
is a (random) finite initial configuration, b,d:W^ x Tq{R*^) —?• [ 0 ; oo) are functions measurable 
with respect to the product cr-algebra =^(M) x ^(ro(M)), and the sequence {..., x_i, xq, xi,...} 
is related to ir]t)te[ 0 -,oo]j as described in Section 1.3.1. We require the processes Ni,N 2 ,r]Q to be 
independent of each other. Equation (31) is understood in the sense that the equality holds a.s. 
for every bounded B G and t > 0. 

As it was said in the preliminaries on Page 6 , we identify a finite configuration with a finite 
simple counting measure, so that a configuration 7 acts as a measure in the following way: 

7(A) = |7nA|, Ag^(M'^). 

We will treat an element of ro(K'^) both as a set and as a counting measure, as long as this 
does not lead to ambiguity. An appearing of a new point will be interpreted as a birth, and a 
disappearing will be interpreted as a death. We will refer to points of r]t as particles. 
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Some authors write d{x,r] \ x) where we write d{x,r]), so that (30) translates to 

LF{rj) = / b{x^r])[F{r]yj x) — F{ri)\dx+ '^^d{x,rj\x){F{rj\x) — F{rj)), (32) 

xeRd 

see e.g. [Pre75], [FKK12b]. 

These settings are formally equivalent: the relation between d and d is given by 

d{x,r]) = d{x,r]\x), ?? G ro(M'^),x G 

or, equivalently, 

d{x, C U x) = d{x, 0, C G ro(M‘^), X G \ e 

The settings used here appeared in [HS78], [GK06], etc. 

We define the cumulative death rate at ^ by 

= (33) 

and the cumulative birth rate by 

BiC)= J b{xX)dx. (34) 

x£R<^ 

2.1 Definition. A (weak) solution of equation (31) is a triple {{r]t)t>o, Ffi, N2), (n,^,P), 
({^t}t>o), where 

(i) (n,#", P) is a probability space, and is an increasing, right-continuous and 

complete filtration of sub - a - algebras of 

(ii) Ni is a Poisson point process on x M_|_ x M+ with intensity dx x ds x du, 

(iii) N2 is a Poisson point process on Z x M+ x M+ with intensity )) x ds x du, 

(iv) r/o is a random ^o’li^sasurable element in ro(M‘^), 

(v) the processes Ni,N 2 and r/o are independent, the processes Ni and N 2 are compatible 
with {^t}t>o, 

(vi) {r]t)t>o is a cadlag ro(K'^)-valued process adapted to ~ d0i 

(vii) all integrals in (31) are well-defined, and 

(viii) equality (31) holds a.s. for all t G [0; 00 ] and all bounded Borel sets B, with {xm}m€Z 
being the sequence related to {'r]t)t>o- 

Note that due to Statement 1.9 item (viii) of this definition is a statement about the joint 
distribution of (t?^), A^i, A^2- 
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Let 


= a{77o, [0; g], C), 7V2(i, [0; (7], C); 

B G G 3S{R+),q G [0;t],i G Z}, 

and let be the completion of 'la^ under P. Note that {^t}t>o is a right-continuous filtration, 
see Remark 1.11. 

2.2 Definition. A solution of (31) is called strong if {r]t)t>o is adapted to > 0). 

2.3 Remark. In the definition above we considered solutions as processes indexed by t G [0; oo). 
The reformulations for the case f G [0;T], 0 < T < oo, are straightforward. This remark applies 
to the results below, too. 

Sometimes only the solution process (that is, {r]t)t>o) will be referred to as a (strong or weak) 
solution, when all the other structures are clear from the context. 

We will say that the existence of strong solution holds, if on any probability space with given 
Ni, N 2 ,riQ, satisfying (i)-(v) of Definition (2.1), there exists a strong solution. 

2.4 Definition. We say that pathwise uniqueness holds for equation (31) and an initial distri¬ 
bution 1 / if, whenever the triples {{r]t)t>o, Ni, N 2 ), (D,^,P), and {{fjt)t>o, Ni, N 2 ), 

(D,^,P), are weak solutions of (31) with P{r]Q = 770 } = 1 and Law{rj) = u, we have 

P{rjt = gtP ^ 1 (that is, the processes 77,77 are indistinguishable). 

We assume that the birth rate h satisfies the following conditions: sublinear growth on the 
second variable in the sense that 


/ b{x,r])dx < Cl It/I -|- C 2 , 


(35) 


and let d satisfy 


We also assume that 


Vttt, G N ; sup d{x,ri) < 00 . 

xSK‘^,|r;|<m 


(36) 


.E|77o| < 00 . (37) 

By a non-random initial condition we understand an initial condition with a distribution, 
concentrated at one point: for some rj' G ro(M'^), P{t 7 o = rj'} = 1. 

From now on, we work on some filtered probability space (D, ({#)}i>o), R). On this 
probability space, the Poisson point processes Ni,N 2 and 770 are defined, so that the whole 
set-up satisfies (i)-(v) of Definition 2.1. 
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Let us now consider the equation 


Vt(.B)= j + ?7o(S), (38) 

Sx (0;t] X [0;oo] 

where h{x,rf) := sup6(a:,^). Note that h satisfies sublinear growth condition (35), if b satisfies 
icv 
it. 

This equation is of the type (31) (with b being the birth rate coefficient, and the zero function 
being the death rate coefficient), and all definitions of existence and uniqueness of solution are 
applicable here. Later a unique solution of (38) will be used as a majorant of a solution to (31). 

2.5 Proposition. Under assumptions (35) and (37), strong existence and pathwise uniqueness 
hold for equation (38). The unique solution {fit)t>o satisfies 

E\fjt\ < oo, t > 0. (39) 

Proof. For co € {f b(x,?]o)dx = 0}, set (t = rjo, an = oo, n ^ N. 

For LO G F := { f b(x, 'rjo)dx > 0}, we define the sequence of random pairs {{an, Co-n)}; where 

(Tn+i = inf{t > 0 : j /[Q.^(^^^^^)](u)diVi(rc,s,u) > 0}+ cr„, fJo = 0, 

X {a„ ;(Tn+d X [0;oo) 

and 

Co — Vo 1 C(7„+1 = Can U {Zn+l{ 

for Zn+1 = {x G ; A^i(x,cJn+i, [0; 6(x, Ccr„)]) > 0}. From (19) it follows that the points Zn are 
uniquely determined almost surely on F. Moreover, > an a.s., and an are hnite a.s. on 

F (particularly because b{x,Can) ^ b{x,r]o)). For uj G F, we define Q = Ca„ for t G [anlcTn+i)- 
Then by induction on n it follows that an is a stopping time for each n G N, and is ^(t„ UiF- 
measurable. By direct substitution we see that {Ct)t>o is a strong solution for (38) on the time 
interval t G [0; lim an)- Although we have not defined what is a solution, or a strong solution, 

n—>-OD 

on a random time interval, we do not discuss it here. Instead we are going to show that 

lim an = oo a.s. (40) 

n—>-cxD 

This relation is evidently true on the complement of F. If P{F) = 0, then (40) is proven. 

If P{F) > 0, define a probability measure on F, Q{A) = , A G ^ Ci F, and define 

= #) n F. 

The process A^i is independent of F, therefore it is a Poisson point process on (F, 5^, Q) 


28 



with the same intensity, compatible with {.y’t}t>o- From now on and until other is specified, we 
work on the filtered probability space {F,^, {^t}t>o, Q)- We use the same symbols for random 
processes and random variables, having in mind that we consider their restrictions to F. 

The process (Ct)te[0; Um cr„) Fas the Markov property, because the process Ni has the strong 

n—>-oo 

Markov property and independent increments. Indeed, conditioning on 


for some xeB} I 


fb{x,CaJdx 

B 


f b(x,Can)dx’ 


thus the chain {C(T„}neZ+ is a Markov chain, and, given {Ca„}neZ+, o'n+i — (Xn are distributed 
exponentially: 

{ -^{o-Ti+i-cr„ >a} I {Ca„}nez+} = exp{-a J b{x,Can)dx}. 

Therefore, the random variables jn = (<7n ~ <^n-i)(/ b{x, Ca^)dx) constitute a sequence of inde- 

pendent random variables exponentially distributed with parameter 1, independent of {Co-„}ne^+- 
Theorem 12.18 in [Kal02] (see Page 22 of this article) implies that (Ci)iG[0; lim an) i® a pure jump 

n—>-oo 

type Markov process. 

The jump rate of (Ci)t€[0; lim < 7 „) is given by 

n—>-oo 

c{a) = / h{x, a)dx. 

Rd 


Condition (35) implies that c(a) < ci|a| + C 2 - Consequently, 

C(C<7„) < Cl|Ca„| +C2 = ClICol +Cin + C2. 

We see that ^ = oo a.s., hence Proposition 12.19 in [Kal02] (given in Section 1.5) 

implies that oo. 

Now, we return again to our initial probability space (12, P”). 

Thus, we have existence of a strong solution. Uniqueness follows by induction on jumps of 
the process. Indeed, let (Ct)f>o be another solution of (38). From (viii) of Definition 2.1 and 
equality 




X (0;c7i) X [0;oo] 

one can see that has a birth before ai} = 0. At the same time, equality 
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x{(Ti }x [0;oo] 

which holds a.s., yields that ( has a birth at the moment ai, and in the same point of space at 
that. Therefore, C coincides with ^ up to ai a.s. Similar reasoning shows that they coincide up 
to an a.s., and, because Un —)• oo a.s., 

P{Ct = Ct for all f > 0} = 1 

Thus, pathwise uniqueness holds. The constructed solution is strong. 

Now we turn our attention to (39). We can write 


101 = \vo\ + - \vo\ > n} 

n=l 

oo 

= \V0\ + ^lWn < t}. 


(41) 


n=l 


Since an = ,■ i we have 


Wn <t} = 


li 


^ fb(x,Cai)dx 


< 


^}c{E 


7* 


^ CilCail +C2 


<n 


n 


7i 


(ci +C2)(|?7 o| +i) 


TT < t} = {Zt - Zo > n}, 


where (Zt) is the Yule process (see Page 19) with birth rate defined as follows: Zt — Zq = n 
when 

n n+1 

7i ^ ^ 7* 


E 


< 


‘<E 


^ (ci + C 2 )(|r?o| +i) ' ^ (ci + C 2 )(|?/o| + i) ’ 
and Zq = \r]Q\. Thus, we have |Ct| < Zt a.s., hence E\Q\ < EZt < oo. □ 

2.6 Theorem. Under assumptions (35)-(37), pathwise uniqueness and strong existenee hold for 
equation (31). The unique solution {r]t) is a pure jump type process satisfying 


E\r]t\ < oo, f > 0. 


(42) 


Proof. Let us define stopping times with respect to > 0}, 0 = 0o E <02 0 03 < •••, 

and the sequence of (random) configurations as follows: as long as 


BiveJ + D{r]0j > 0 , 
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we set 


9n+l — ^n+l ^n+1 “I" 

6*^+1 = inf{t > 0 : J /[ 0 ;b( 3 ,,,,e^)](«)diVi(x, s, u) > 0 }, 

R'* X (dn ]dn+t] X [0;oo) 

On+1 = mf{t > 0 ; j I{^.(,^g^}I[0.dixurigJ]iv)dN2{i,r,v) > 0 }, 

(dn-,dn+t]x[0-,Oo) 

%„+i = Ve„ U {Zn+i} if (9^+1 < 9i+i, where {zn+i} = {z G M"* : Ni{z,9n + (9^+i,M+) > 0}; 
^?e„+i = r?e„ \ {zn+i} if 6'n+i > 6'n+i> where {zn+i} = {xi G r]e„ : N 2 {i, 9n + 6'^+i,M+) > 0 }; the 
configuration = r/o is the initial condition of (31), r]t = r/e„ for t G [6*„;6*„+i), {xi} is the 
sequence related to {rjt)t>o- Note that 

P{9t+i = 6*^+1 for some n \ B{rjeJ + D{rjgJ > 0} = 0, 

the points Zn are a.s. uniquely determined, and 

P{zn+i G I < 0^+1} = 0. 


If for some n 

Bimn) + Pimn) = 0 , 

then we set 9n+k = co, /c G N, and rjt = rje^, t > 9n- 

As in the proof of Proposition 2.5, [rjt) is a strong solution of (31), t G [O;hm„0„). 

Random variables 9n,n G N, are stopping times with respect to the filtration {^tP P 0}. 
Using the strong Markov property of a Poisson point process, we see that, on {0„ < oo}, the 
conditional distribution of given is exp( f b(x, r]g^)dx), and the conditional distribution 

of given ^g^ is exp( ^ d{x,'r]g^)). In particular, 9^,9'^ > 0, n G N, and the process {rjt) 
is of pure jump type. 

Similarly to the proof of Proposition 2.5, one can show by induction on n that equation (31) 
has a unique solution on Namely, each two solutions coincide on [0; 0,^] a.s. Thus, any 

solution coincides with {rjt) a.s. for all t G [O;0„]. 

Now we will show that 9n ^ oo a.s. as n —)• oo. Denote by the moment of the k-th. birth. 
It is sufficient to show that ^ oo, /c ^ oo, because only finitely many deaths may occur 
between any two births, since there are only finitely particles. By induction on k' one may see 
that C {cTjliep}, where at are the moments of births of ijjt)t>o^ the solution of (38), 

and rjt C rjt for all t G [O;hmjj0„). For instance, let us show that (%)t>o has a birth at 9[. We 
have rigi^_ D tjq = rjQ, and r]gi^_ C rjt |t=o= Vo, hence for all x G 
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b{x,r]g,^_) > > 6(x,r/0/_) 

The latter implies that at time moment 6[ a birth occurs for the process (%)t>o in the same 
point. Hence, C rjg/, and we can go on. Since cifc —)• oo as A: ^ oo, we also have 0^ —)■ oo, 
and therefore ^ oo, n —)■ oo. 

Since rjt C rj^ a.s.. Proposition 2.5 implies (42). □ 

In particular, for any time t the integral 

R‘*x(0;t]x[0;oo] 

is finite a.s. 

2.7 Remark. Let r/o be a non-random initial condition, rjo = a, a ^ ro(K'^). The solution of 
(31) with r]Q = a will be denoted as (r]{a,t))t>o- Let be the push-forward of P under the 
mapping 


n 9 w i-^-( 77 ( 0 ;, •)) G T]. (43) 

From the proof one may derive that, for hxed cu G fl, constructed unique solution is jointly 
measurable in {t,a). Thus, the family {Pa} of probability measures on T] is measur¬ 
able in a. We will often use formulations related to the probability space (L^rolK^*) [Oj T]),Pa) 

in this case, coordinate mappings will be denoted by rjt, 

r]t{x) = x{t), X G Zlpg(Kd)[0;T]. 

The processes {'i]t)te[ 0 ;T] {r]{a, •))te[ 0 ;r] have the same law (under Pa and P, respectively). 

As one would expect, the family of measures {Pa, ct £ ro(M'^)} is a Markov process, or a Markov 
family of probability measures; see Theorem 2.15 below. For a measure on ro(K'^), we define 

-P/i = y PaKdct). 

We denote by the expectation under 

2.8 Remark. Let bi, di be another pair of birth and death coefficients, satisfying all conditions 
imposed on b and d. Consider a unique solution (fjt) of (31) with coefficients 61 , di instead 
of b, d, but with the same initial condition r]Q and all the other underlying structures. If for 
all ( £ D, where D G .^(ro(M'^)) , 6 i(-,C) = b{-X), '^i(gC) = d{-X), then fjt = rjt for all 
t < inf{s > 0 '■ Tjs ^ D} = inf{s > 0 ■ fjs ^ D}. This may be proven in the same way as the 
theorem above. 
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2.9 Remark. Assume that all the conditions of Theorem 2.6 are fulfilled except Condition 
(37). Then we could not claim that (42) holds. However, other conclusions of the theorem 
would hold. We are mostly interested in the case of a non-random initial condition, therefore 
we do not discuss the case when (42) is not satisfied. 

2.10 Remark. We solved equation (31) w-wisely. As a consequence, there is a functional 
dependence of the solution process and the “input”: the process {Tjt)t>o is some function of r]Q, 
Ni and A^ 2 - Note that 9n and Zn from the proof of Theorem 2.6 are measurable functions of 
rjQ, Ni and N 2 in the sense that, e.g., 0i = Fi{rio, Ni, N 2 ) a.s. for a measurable Fi : ro(M'^) x 
r(M‘^ X M+ X M+) X r(Z"' X M+ X M+) ^ M+. 

2.11 Proposition. If {r}t)t>o is a solution to equation (31), then the inequality 

E\r]t\ < (c2t-h.E|r/o|)e‘'i* 


holds for all t > 0. 

Proof. We already know that E\r]t\ is finite. Since r]t satisfies equation (31) we have 


m{B)= j IioXx,ris-)]iu)dNi{x,s,u) 

Sx (0;t] X [0;oo] 


Zx (0;t] X [0;oo) 

j I[ 0 M^,ris-)]i'>^)dNi{x,s,u)+r]o{B). 

Bx (0;£] X [0;oo] 

For B = taking expectation in the last inequality, we obtain 

E\r]t\ = <E j I[Q.h(^^^^_^_)]{u)dNi{x,s,u) + Er]o{R^) = 

R^x(0;£]x[0;oo] 

= E J I[o.^^^^_^'^{u)dxdsdu + Er]Q{R'^) = E J b{x,ris-)dxds + EriQ{R'^). 

IR'^x(0;t]x[0;oo] R‘*x(0;t] 

Since 77 is a solution of (31), we have for all s G [0;f] almost surely = 'r]s- Consequently, 
E\r]s-\ = E\ris\. Applying this and (35), we see that 

Er]t{R‘^) < E J {ci\r]s-\ + C 2 )ds + Erjoi'R'^) = Cl J E\r]s\ds + C 2 t + ErjofR'^), 

(0;t] (0;t] 

so the statement of the lemma follows from (37) and Gronwall’s inequality. □ 
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2.12 Definition. We say that joint uniqueness in law holds for equation (31) with an initial 
distribution n if any two (weak) solutions ({rit), Ni, N 2 ) and ((rjt)', N[, N 2 ) of (31), Law{rjo) = 
Law{{riQ)') = V, have the same joint distribution: 


Law{{r,t),N^,N 2 ) = Law{{7^t)' ,N[,N! 2 ). 

The following corollary is a consequence of Theorem 2.6 and Remark 2.10 . 

2.13 Corollary. Joint uniqueness in law holds for equation (31) with initial distribution v 
satisfying 


/ \^\i'{d'y) < 00 . 

2.14 Remark. We note here that altering the order of the initial configuration does not change 
the law of the solution. We could replace the lexicographical order with any other. To see this, 
note that if ? is a permutation of Z (that is, <;■ ; Z —)• Z is a bijection), then the process N 2 
defined by 

N2{K,R,V) = N2{^K,R,V), Kcz,i?,l/G^(M+), (44) 

has the same law as iV 2 , and is adapted to too. Therefore, solutions of (31) and of 

(31) with N 2 being replaced by N 2 have the same law. But replacing N 2 with N 2 in equation 
(31) is equivalent to replacing ...,Xo,xi, ...} with 

Let be a distribution on ro(K'^), and let T > 0. Denote by ^{v,b,d,T) the law of the 
restriction (??t)tg[o;r] of the unique solution {'r]t)t>o to (31) with an initial condition distributed 
according to n. Note that ^{n,b,d,T) is a distribution on Dp^(gd)([0; T]). As usually, the 
Markov property of a solution follows from uniqueness. 

2.15 Theorem. The unique solution {'r]t)t€[ 0 ;T] of (31) ® Markov proeess. 

Proof. Take arbitrary t' < t, t',t G [0;T]. Consider the equation 


^tiB)= j /[o;b(,i,,^^_)](u)(iiVi(x,s,u) 

X [0;oo] 

A v) + ht' {B), 


Zx X [0;oo) 


(45) 


where the sequence {x(} is related to the process (^s)sg[o;t]) = hs- The unique solution of 

(45) is As in the proof of Theorem 2.6 we can see that {'Hs)s&[t';t] is measurable with 
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respect to the filtration generated by the random variables Ni{B, [s; g], U), N 2 {i, [s; q], U), and 
where B G iGZ, t'<s<q<t, UG Poisson point process have 

independent increments, hence 

P{{Vt)te[s;T] £ U I ^s} = P{{Vt)te[s-,T] £ U I Vs} 

almost surely. Furthermore, using arguments similar to those in Remark 2.14, we can conclude 
that iris)s£[t'-,t] is distributed according to — t'), where i^t' is the distribution of r/^'. 

□ 

The following theorem sums up the results we have obtained so far. 

2.16 Theorem. Under assumptions (35), (36), (37), equation (31) has a unique solution. This 
solution is a pure jump type Markov process. The family of push-forward measures {Pa, a: G 
ro(K'^)} defined in Remark 2.1 forms a Markov process, or a Markov family of probability mea¬ 
sures, on Zlp^(]Rd)[0;oo). 

Proof. The statement is a consequence of Theorem 2.6, Remark 2.7 and Theorem 2.15. In 
particular, the Markov property of {Pa,a G ro(K'^)} follows from the statement given in the 
last sentence of the proof of Theorem 2.15. □ 

We call the unique solution of (31) (or, sometimes, the corresponding family of measures on 
UpofR'*) [Oj oo)) ^ (spatial) birth-and-death Markov process. 

2.17 Remark. We note that d does not need to be dehned on the whole space x ro(M'^). 
The equation makes sense even if d{x,r]) is dehned on {{x,r]) \ x G rj}. Of course, any such 
function may be extended to a function on x ro(K‘^). 

2.1 Continuous dependence on initial conditions 

In order to prove the continuity of the distribution of the solution of (31) with respect to initial 
conditions, we make the following continuity assumptions on b and d. 

2.18 Continuity assnmptions. Let b,d be continuous with respect to both arguments. Fur¬ 
thermore, let the map 

ro(M‘^) 9 7?^ 6(-,r/) G L\R^) 

be continuous. 

In light of Remark 2.17, let us explain what we understand by continuity of d when d{x,ri) 
is dehned only on {{x, p) \ x G r]}. We require that, whenever Pn ^ V and 9 Zjj —)• x G r/, we 
also have d{zn,r]n) d{x,r]). Similar condition appeared in [HS78, Theorem 3.1]. 
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2.19 Theorem. Let the birth and death coefficients b and d satisfy the above continuity assump¬ 
tions 2.18. Then for every T > 0 the map 


ro(M'^) 9 a !-)■ Law{ri{a, ■).,■£ (0;T]}, 

which assigns to a non-random initial condition r]Q = a the law of the solution of equation (31) 
stopped at time T, is continuous. 

Remark. We mean continuity in in the space of measures on Page 13. 

Proof. Denote by 77 ( 0 , •) the solution of (31), started from a. Let —?• a, an, a € ro(K'^), 
a = {xo,X-i ,..., xo ^ x_i ^ ... ^ x_|q|+i. With no loss in generality we assume that 

\an\ = |ck|) u G N. By Lemma 1.5 we can label elements of = {xl^\x^ffl, •••, so 

(n) 

that x_f X-i, 7 = 0,..., |q;| — 1. Taking into account Remark 2.14, we can assume 


M) ^ 


^0 




^ x 


(n) 

-|a|+l 


(46) 


without loss of generality (in the sense that we do not have to use lexicographical order; not in 
the sense that we can make x^^\x^ffl,... satisfy (46) with the lexicographical order). 

We will show that 


sup dist{r]{a,t),r]{an,t)) A-0, n —)• 00 . (47) 

te[ 0 ;r] 

Let be the moments of jumps of process 77 ( 0 , •). Without loss of generality, assume 

that d{x,a) > 0, x G a, and ||ft(-,a)||Li > 0, := L^(M'^) (if some of these inequalities are not 

fulfilled, the following reasonings should be changed insignificantly). 

Depending on whether a birth or a death occurs at 61 , we have either 

W({a;i} X { 01 } X [0;6(xi,r7o)]) = 1 (48) 


or for some X-k G a 

N 2 {{-k} X {0i} X [0;(i(x_fc,a)]) = 1. 

The probability of last two equalities holding simultaneously is zero, hence we can neglect this 
event. In both cases W(a;i, {0i}, { 6 (xi, a)}) = 0, N2{—k,{9i},{d{x-k,a)}) = 0 a.s. We also 
have 

X [0; 9i) X [0; b{x, a)]) = 0, 

and for all j G 0,1,..., |q;| — 1 

N 2 {{-j} X [O;0i) X [0;d{x-j,a)]) =0. 
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Denote 

m := b{xi, a) A min{(i(x, a) : x G a} A ||6(-, a)||Li A 1 

and fix e > 0. Let (5i > 0 be so small that for u G ro(K'^), ly = {xq, x'_;^,x'_|^|_|_j^}, \x-j—x'_j\ < 
(5i the inequalities 

\d{x'_j, u) — d{x-j, a)I < em, 

and 

\\b{-,iy) - b{-,a)\\Li < em 

hold. Then we may estimate 

P{ j > 1} < e. 

K'ix[0;6»i)x[0;oo] 

and 

P{ j I{xmeiy}Ilo-,d{xm,u)]iv)dN 2 ii,r,v) >1} <E\a 

Zx[O;0i)x[O;oo] 

Indeed, the random variable 


(49) 

(50) 
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X [0;^) X [0;oo] 


'^[0;0V{6(a:,L')—6(x,a})] (^5 


(51) 


is exponentially distributed with parameter ||(6(-,z^) — 6(-,a))+||^i < e||6 (-, q;)| l^i. By Lemma 


1.13, 


P{e < Oi} < 


^\\b{-,a)\\p 

IIK->a)llLi 


= e, 


(52) 


which is exactly (49). Likewise, (50) follows. 

Similarly, the probability that the same event as for r]{a, •) occurs at time 9i for ri{u, ■) is 
high. Indeed, assume, for example, that a birth occurs at 9i, that is to say that (48) holds. 
Once more using Lemma 1.13 we get 


F{JV,({x,) X {«i) X [O;(,(j,i,01) =0} < 

The case of death occurring at 6 i may be analyzed in the same way. 
From inequalities (9) and (10) we may deduce 


sup dist{T]{a,t),r]{an,t)) ^ 0,n ^ oo. (53) 

te( 0 ;ei] 
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Proceeding in the same manner we may extend this to 


sup dist{r]{a,t),r]{an,t)) A-0,n ^ oo, (54) 

te(0;6»„] 

particularly because of the strong Markov property of a Poisson point process. In fact, with 
high probability the processes •) and r]{a, ■) change up to time 9^ in the same way in the 
following sense: births occur in the same places at the same time moments. Deaths occur at the 
same time moments, and when a point is deleted from r]{a, ■), then its counterpart is deleted 
from r]{an, ■ )■ 

Since 0^ —)• oo, we get (47). □ 

2.20 Remark. In fact, we have proved an even stronger statement. Namely, take — ?• a. 

Then there exist processes (^t”^)tG[ 0 ;T] such that {Ci"^^)tG[o-,T] = (^(®n, ^))iG[ 0 ;T] 

sup A- 0, n^oo. 

te[o-,T] 

Thus, Law{r]{a,-),- £ (0;T]} and Law{r]{an,-),- £ (0;T]} are close in the space of measures 
over Dyq , even when is considered as topological space equipped with the uniform topology 
(induced by metric dist), and not with the Skorokhod topology. 

2.2 The martingale problem 

Now we briefly discuss the martingale problem associated with L defined in (30). Let C'b(ro(K'^)) 
be the space of all bounded continuous functions on Po(K'^). We equip C'f,(ro(K'^)) with the 
supremum norm. 

2.21 Definition. A probability measure Q on (ZlrQ[0; oo), .^(Urg [0; oo))) is called a solution 
to the local martingale problem associated with L if 

t 

m/ = f{y{t)) - fiy{0)) - J Lf{y{s-))ds, 0<t <oo, 

0 

is a local martingale for every / £ C'b(ro). Here y is the coordinate mapping, y{t){u!) = 
u £ Dro[0;oo), is the completion of a{y{s),0 < s <t) under Q. 

Thus, we require to be a local martingale under Q with respect to {^t}t>o- Note that 
L can be considered as a bounded operator on C'b(ro(K'^)). 
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2.22 Proposition. Let {r]{a,t))t>o be a solution to (31). Then for every f G C'(ro) the process 

t 

= fivia, t)) - fir]{a, t)) - j Lf{r]{a, s-))ds (55) 

0 

is a local martingale under P with respect to {^t}t>o- 

Proof. In this proof Q will stand for 'q{a,t). Denote = inf{f > 0 ; |Ct| > n or ft ^ 
[—n]nY‘}. Clearly, r^, n G N, is a stopping time and Tn ^ oo a.s. Let ff = ftArn- We want to 
show that m/) t>o is a martingale, where 

t 

= f{Q) - /(CD - I Lf{Q_)ds. (56) 

0 

The process {ft)t>o satisfies 


Ci= J] [Cs-Cs-] + Co. (57) 

In the above equality as well as in few other places throughout this proof we treat elements of 
ro(K'^) as measures rather than as configurations. Since (ft) is of the pure jump type, the sum 
on the right-hand side of (57) is a.s. finite. Consequently we have 

/(CD - /(CoD = E - /(CD)] 

= j ifiCs) - fiCs-)]I{s<T„}I[ 0 Mx,Cs-)]i'^)dNi{x,s,u) (58) 

Sx (0;t] X [0;oo] 

J h^i&(s-}[fiCs) - f{Cs-)]I{s<T„}I[o-,di^i,(:s-)]iv)dN2{i,s,v). 

Zx (0;t] X [0;cxd] 

Note that Cs = Cs- U x a.s. in the first summand on the right-hand side of (58), and 
Cs = Cs- \ Xi a.s. in the second summand. Now, we may write 

t 

J /{,<,„}L/(Cs)ds = 

0 

hs<r„}I[o-Mx,Cs-)]i^)[fi(s- u x) - f{Cs-)]dxduds- (59) 

0 xGR^,u>0 
t 

J j I{s<T„}I[o-,d{xXs-))]iv)[fiCs- \x)- f{Cs-)]Cs-{dx)dvds. 

0 xGM^,u>0 
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Functions b,d{-,-) and / are bounded on x {a : |q;| < n and a C [—n;n]'^} and {a : 
|a| < n and a C [—n;n]‘^} respectively by a constant C > 0. Now, for a predictable bounded 
processes (7s(x, u))o<s<t and {/3s{x,v))Q<s<ti the processes 


J I{s<Tn}ls{x,u)[dNi{x,s,u) - dxdsdu], 

Bx(0;t]x[0;C] 

- #{di)dsdv]. 


Zx{0;t]x[0-,C] 

are martingales. Observe that 


J I{s<Tn}h^i(^(:s-}Mxi,v)i^{di)dsdv = J I[s<Tn}Ps{x,v)Cs-{dx)dsdv 

Zx(0;t]x[0;C] Zx(0;t] x [0;C] 

Taking 


-fs{x,u) = I[0;b(a;,C^_)](u)[/(Cs- Ux) - /(C_)], 

(ds{x,v)) = I[o.d(x,(;,-)]{v)[f{Cs-\x) - fiCs-)], 

we see that the difference on the right hand side of (56) is a martingale because of (58) and 
(59). □ 

2.23 Corollary. The unique solution of (31) induces a solution of the martingale problem 2.21. 

2.24 Remark. Since y{s) = y{s—) Pa - a.s., the process 

t 

fivit)) - /(y(0)) - J Lf{y{s))ds, 0<t <oo, 

0 


is a local martingale, too. 


2.3 Birth rate without sublinear growth condition 

In this section we will consider equation (31) with the a birth rate coefficient that does not 
satisfy the sublinear growth condition (35). 

Instead, we assume only that 


sup b{x, y) < oo. 


(60) 
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Under this assumption we can not guarantee existence of solution on the whole line [0; oo) 
or even on a finite interval [0; T], It is possible that infinitely many points appear in finite time. 

We would like to show that a unique solution exists up to an explosion time, maybe finite. 
Consider birth and death coefficients 


bn{x,r]) = 

dn{x,r]) = d{x,ri)I^\^\<n}- 


(61) 


Functions bn,dn are bounded, so equation (31) with birth rate coefficient bn and death rate 
coefficient dn has a unique solution by Theorem 2.6. Remark 2.8 provides the existence and 
uniqueness of solution to (31) (with birth and death rate coefficients b and d, respectively) up 
to the (random stopping) time = inf{s > 0 ; \rjs\ > n}. Clearly, r^+i > if —)• oo a.s., 
then we have existence and uniqueness for (31); if t < oo with positive probability, then 
we have an explosion. However, existence and uniqueness hold up to explosion time r. When 
we have an explosion we say that the solution blows up. 


2.4 Coupling 

Here we discuss the coupling of two birth-and-death processes. The theorem we prove here will 
be used in the sequel. As a matter of fact, we have already used the coupling technique in the 
proof of Theorem 2.6. 

Consider two equations of the form (31), 


(H) = 




Bx{0;t] X [0;oo] 


(H), fc = 1,2, 


(62) 


Zx (0;t] X [0;oo) 


where t G [0;r] and is the sequence related to (^i^^)te[ 0 ;r]- 

(k) 

Assume that initial conditions Q ^ and coefficients b^, dfc satisfy the conditions of Theorem 
2.6. Let (^t^^)t 6 [ 0 ;r] b® the unique strong solutions. 

2.25 Th 

1 9 

r? C ry , 


2.25 Theorem. Assume that almost surely C and for any two finite configurations 


and 


bi{x,r]^) < b 2 {x,r]‘^), x G 

di{x,r]^) > d 2 ix,r]‘^), x G ry^. 


(63) 
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Then there exists a cadlag Fq 
have the same law and 


-valued proeess {vt)tG[o-,T] such that {r]t)te[o-,T] and ^)te[o;r] 


ritCCi^\ tG[0;T]. 


(64) 


Proof. Let {..., re®, ...} be the sequence related to 


( 2 )n 


iG[0;T]- 


equation 


Consider the 


VtiB) = 


Ilo-Mx,,j,_)]{u)dNi{x, s,u) 


Bx (0;£] X [0;oo] 




(65) 


Zx (0;t] X [0;oo) 


Note that here {x) '} is related to (^) JtG[ 0 ;T] not to {'i]t)t£[o-,T]- Thus (65) is not an equation 
of form (31). Nonetheless, the existence of a unique solution can be shown in the same way as 
in the proof of Theorem 2.6. Denote the unique strong solution of (65) by {r]t)t&[ 0 ;T]- 

(2i') 

Denote by {r^lmeN the moments of jumps of {'nt)te[o-,T] and {Q )iG[ 0 ;T], 0 < n < r 2 < rs < 

( 2 ) 

.... More precisely, a time t G {rm}mGN iff at least one of the processes {f]t)t£[ 0 ;T] and {Q Ote[ 0 ;r] 
jumps at time t. 

We will show by induction that each moment of birth for ir]t)te[o-,T] is a moment of birth for 

f2) i‘2') 

{Q Jte[ 0 ;r] too, and each moment of death for JtG[ 0 ;T] is a moment of death for (??i)tg[o;r] if 

the dying point is in (?7i)iG[0;T]- Moreover, in both cases the birth or the death occurs at exactly 

the same point. Here a moment of birth is a random time at which a new point appears, a 

moment of death is a random time at which a point disappears from the configuration. The 

statement formulated above is in fact equivalent to (64). 

Here we deal only with the base case, the induction step is done in the same way. We have 

nothing to show if ri is a moment of a birth of (^) Jte[0;r] or a moment of death of {r}t)t&[o-,T]- 

Assume that a new point is born for (??t)te[o;T] D) 


Vri \ rjTi- = {a:i}. 

The process {'i]t)te[o-,T] satisfies (65), therefore {d}, [0; 6 fc(xi, ??t-i-)]) = 1- Since 

„ _ g(1) ^ g( 2) _ f(2) 

'm-— SO ^ so — Sti-) 


by (63) 


[ 0 ; 6 fc(a:i,?®_)]) = 1 , 


42 



hence 

= {xi}. 

(2i') 

The case when T 2 is a moment of death for {Q Jte[ 0 ;r] is analyzed analogously. 

It remains to show that {rjt)te[ 0 ;T] ^i^cl (^t^^)t 6 [ 0 ;r] have the same law. We mentioned above 
that formally equation (65) is not of the form (31), so we can not directly apply the unique- 
ness in law result. However, since r]t G Q a.s., t G [0;T], we can still consider (65) as 
an equation of the form (31). Indeed, let {..., y_i, yo, yi,...} be the sequence related to ijt- 
We have {y_|^(i)|_^^,..., y_i, yo, 2 / 1 , •••} C {x_|^( 2 )|_^^,..., xf\ ...}. There exists an in¬ 
jection ? : { —+ 1,..., 0,1,...} + 1,..., 0,1,...} such that y^(j) = Xj. Denote 

9i = inf{s > 0 : y* G r]s}. Note that 9i is a stopping time with respect to {^t}- Define a Poisson 
point process N 2 by 

/V 2 ({f} xRxV) = iV 2 ({f} xRxV), i£Z,Rc[0;9i],V C M+, 


and 

Af2({f} xRxV) = iV2({?(f)} xRxV), i£Z,Rc ( 0 *; 00), V C M+. 

The process N 2 is {^t}-adapted. One can see that ('f]t)te[ 0 ;T] is the unique solution of equation 
(31) with N 2 replaced by N 2 . Hence (7]t)te[0;T] = 

2.5 Related semigroup of operators 

We say now a few words about the semigroup of operators related to the unique solution of (31). 
We write ri{a,t) for a unique solution of (31), started from a G ro(K'^). We want to define an 
operator St by 

Stf{a) = Ef{r]{a,t)) {= Eafivit))) (66) 

for an appropriate class of functions. Unfortunately, it seems difficult to make St a Co-semigroup 
on some functional Banach space for general b,d satisfying (35) and (36). 

We start with the case when the cumulative birth and death rates are bounded. Let Cb = 
C;,(ro(K'^)) be the space of all bounded continuous functions on ro(K'^). It becomes a Banach 
space once it is equipped with the supremum norm. We assume the existence of a constant 
C > 0 such that for all () G ro(K'^) 


|H(C)| + |D(C)| <C, (67) 

where B and D are defined in (33) and (34). Formula (30) defines then a bounded operator 
L : Cb ^ Cb, and we will show that St coincides with For / G Cb, the function Stf is 
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bounded and continuous. Boundedness is a consequence of the boundedness of /, and continuity 
of Stf follows from Remark 2.20. Indeed, let an a, = r]{an,t) and 

dist{ri{a,t),^j!^^) ^ 0, n ^ oo. 

Unlike r(M'^), the space ro(M'^) is a fi-compact space. Consequently, for all e > 0 there exists a 
compact Ke C ro(K'^) such that for large enough n 

Also, for fixed (5 > 0 and for large enough n 

P{dist{ri{a,t),^^'^) < (5} > 1 — A 

Fix e > 0. There exists 5e G (0;e) such that |/(/3) — f{'y)\ < e whenever dist{f3,'y) < 6^, 
/3 ,7 G Kg. We have for large enough n 

\E[f(r)(a,t)) - /(Cj”^)]| 

< Elf(r}(a,t)) - /(Ci^”^)|/{r/(a,f) G Kg, G Kg,dist(?](a,t),^^'"^) < (5£} 

+2((5£ + e)||/|| < e + 2((5£ + e)||/||, 

where ||/|| = sup^gp^^l^d) 1/(01- Letting e ^ 0, we see that 

Ef(Tj(an,t)) = ^ Ef{rj{a,t)). 

Thus, Stf is continuous (note that the continuity of Stf does not follow from Theorem 2.19 
alone, since for a fixed t G [0;T] the functional ^ x{t) G M is not continuous 

in the Skorokhod topology). Furthermore, since for small t and for all A G 


and for x £ a 


P{ri{a, t) = a} = 1 — t[B{a) + D{a)] + o{t), 

P{ri{a, f) = a U {y} for some y £ A} = t / b{y, a)dy + o{t), 

y&A 

P{ri{a, t) = a \ {x}} = td{x, a) + o{t), 


( 68 ) 

(69) 

(70) 


we may estimate 
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\Stf{a) - f{a)\ < t[B{a) + D{a)] \ \f\ \ +o(t)||/|| < C\\f\\t + o{t). 


Therefore, (66) defines a Cq semigroup on Cb- Its generator 


Lf{a) 


lim 
t — 


Stfja) 

t 


lim 
t —^0+ 


/ 


b{x, a)[f{a U x) — f{a)]dx + d{x, a){f{a \x) — /(a)) + o{t) 

xGa 


Thus, St = and we have proved the following 


= Lf{a). 


2.1 Proposition. Assume that (67) is fulfilled. Then the family of operators {St,t > 0) on Cb 
defined in (66) constitutes a Co-semigroup. Its generator coincides with L given in (30). 

Now we turn out attention to general b,d satisfying (35) and (36) but not necessarily (67). 
The family of operators {St)t>o still constitutes a semigroup, however it does not have to be 
strongly continuous anymore. Consider truncated birth and death coefficients (61) and corre¬ 
sponding process ri'^{a,t). Remark 2.8 implies that r]'^{a,t) = rj(a,t) for all t G [0;rn], where 


Tn = inf{s > 0 : \r]{a, s)| > n}. (71) 

Growth condition (35) implies that —)• oo for any a G ro(M'^). 

Truncated coefficients bn,dn satisfy (67) and 

si^'>f{a) = Ef{r^^^\a,t)) (72) 

defines a Cq - semigroup on Cb. In particular, for all a G ro(M'^) 

L(»)/(„)= lin, En„l-Ha,t)) - f(o) 

t^o-e t 

where is operator defined as in (30) but with bn, dn instead of b, d. Letting n —)• oo we get, 
for fixed a and /, 


Lf{a) 


lim 
t —^0-|- 


Ef{r]{afi)) - f{a) 
t 


lim 
t — 


Stfja) - /(«) 

t 


(73) 


Taking limit by n is possible: for n > |a| +2, r]^'^\a,t) satisfies (68), (69) and (70), therefore 
rija,t) satisfies (68), (69) and (70), too. Thus, we have 


2.2 Proposition. Let b and d satisfy (35) and (36) but not necessarily (67). Then the family of 
operators jSt, t >0) constitutes a semigroup on Cb which does not have to he strongly continuous. 
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However, for every a G ro(K'^) and f £ Cb we have (73). 

Formula (73) gives us the formal relation of {r]{a,t))t>o to the operator L. Of course, for 
fixed / the convergence in (73) does not have to be uniform in a. 

2.3 Remark. The question about the construction of a semigroup acting on some class of 
probability measures on ro(M'^) is yet to be studied. 
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